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Abstract 

In this paper the difference in the asymptotic dynamics between the nonlocal and 
local two-dimensional Swift-Hohenberg models is investigated. It is shown that the 
bounds for the dimensions of the global attractors for the nonlocal and local Swift- 
Hohenberg models differ by an absolute constant, which depends only on the Rayleigh 
number, and upper and lower bounds of the kernel of the nonlocal nonlinearity. Even 
when this kernel of the nonlocal operator is a constant function, the dimension bounds 
of the global attractors still differ by an absolute constant depending on the Rayleigh 
number. 
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1 Introduction 



Fluid convection due to density gradients arises in geophysical fluid flows in the atmo- 
sphere, oceans and the earth's mantle. The Rayleigh-Benard convection is a prototypical 
model for fluid convection, aiming at predicting spatio-temporal convection patterns. The 
mathematical model for the Rayleigh-Benard convection involves nonlinear Navier-Stokes 
partial differential equations coupled with the temperature equation. When the Rayleigh 
number is near the onset of the convection, the Rayleigh-Benard convection model may be 
approximately reduced to an amplitude or order parameter equation, as derived by Swift 
and Hohenberg (p5||). 

In the current literature, most work on the Swift-Hohenberg model deals with the 
following one-dimensional equation for w{x,t), which is a localized, one-dimensionalized 
version of the model originally derived by Swift and Hohenberg ( [|15| ) , 

Wt = nw-{l + dxxfw-w^. (1) 

The cubic term is used as an approximation of a nonlocal integral term. For the (local) 
one-dimensional Swift-Hohenberg equation (|^), there has been some recent research on 
propagating or steady patterns (e.g., |^], [^). Mielke and Schneider(||l^) proved the 
existence of the global attractor in a weighted Sobolev space on the whole real line. Hsieh 
et al. ([^, [^) remarked that the elemental instability mechanism is the negative diffusion 
term —Wxx- 

Roberts ([^], [13|) recently re-examined the rationale for using the Swift-Hohenberg 



model as a reliable model of the spatial pattern evolution in specific physical systems. He 
argued that, although the localization approximation used in (|l]) makes some sense in the 
one-dimensional case, this approximation is deficient in the two-dimensional convection 
problem and one should use the nonlocal Swift-Hohenberg model (p|, ||], H): 



ut 



liu-{l + Afu -uj^ {x - + {y - iiY)u\i, r,, t)dCdv, (2) 



where u = u{x, y, t) is the unknown amplitude function, ^ measures the difference of the 
Rayleigh number from its critical onset value, A = dxx + dyy is the Laplace operator, and 
G(r) is a given radially symmetric function (r = + y^). The equation is defined for 
t > and {x,y) S D, where D is a bounded planar domain with smooth boundary dD. 
The two-dimensional version of the local Swift-Hohenberg equation for u{x, y, t) is 

Ut = /iu - (1 + A)^M - u^. (3) 

,3 ; 



Here u is used to approximate the nonlocal term in (y). 



Roberts ([0], [13|) noted that the range of Fourier harmonics generated by the non- 
linearities is fundamentally different in two-dimensions than in one-dimension. This dif- 
ference requires a more sophisticated treatment of two-dimensional convection problem, 
which leads to nonlocal nonlinearity in the Swift-Hohenberg model. He also argued that 
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nonlocal operators naturally appear in systematic derivation of simplified models for pat- 
tern evolution, and nonlocal operators also permit symmetries which are consisitent with 
physical considerations. 

In this paper, we discuss the difference between nonlocal and local two-dimensional 
Swift-Hohenberg models (^), (^), from a viewpoint of asymptotic dynamics. We show that 
the bounds for the dimensions of the global attractors for the nonlocal and local Swift- 
Hohenberg models differ by an absolute constant, which depends only on the the Rayleigh 
number, and upper and lower bounds of the kernel of the nonlocal nonlinearity. Even 
when this kernel is a constant function, the dimension bounds of the global attractors still 
differ by a constant depending on the Rayleigh number. In §2 and §3, we will consider the 
nonlocal and local Swift-Hohenberg models, respectively. Finally in §4, we summarize the 
results. 



2 Nonlocal Swift-Hohenberg Model 

In this section, we discuss the global attractor and its dimension estimate for the nonlocal 
Swift-Hohenberg model (^). In the following we use the abbreviations = L'^{D), L°° = 
L°°{D), = H^{D) and Hq = Hq{D) {k is a non-negative integer) for the standard 
Sobolev spaces. Let (•,•), || • || = || • ||2 denote the standard inner product and norm in 
L^, respectively. The norm for Hq is || • \\fjk. Due to the Poincare inequality, ||L)'^n|| is an 

equivalent norm in Hq. 

We rewrite the two-dimensional nonlocal Swift-Hohenberg equation (|2|) as 

ut + au + 2Au + A\ + n {x - 0^ + iv - V, t)d^dr] = 0, (4) 

where a = 1 — ii. This equation is supplemented with the initial condition 

u{x,y,0) =uo{x,y), (5) 

and the boundary conditions 

du . . 

u\aD = 0, ^|9D = 0, (6) 

where n denotes the unit outward normal vector of the boundary dD. 

In this paper, we assume the following conditions for every t > and {x, y) G D, 

0<b< G{^x^ + y2) < a, and G, VG, AG G L'^iD), (7) 

where a,b > are some positive constants and V = {dx,dy) is the gradient operator. 
Denote Ki = ||VG||oo and K2 = ||AG||oo- 

To study the global attractor, we need to derive some a priori estimates about solutions. 
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Lemma 1 Suppose u is a solution of Then u is uniformly (in time) bounded, 

and the following estimates hold for t > 



and thus 



where R = ^/^. 



\u{x,y,t)f < \\uo{x,y)fexp{-2fit) + ^, (8) 



lim supt^^^\\u{x, y,t)\\ < ^ = R, (9) 



Proof. Taking the inner product of (Q) with u, we have 

^^||uf + \\Auf + 2{Au,u) + a\\uf 



Note that 



+(n2, G{^{x-Cy + {y-VyW{t ri)didi^) = 0. (10) 
2\{Au,u)\ < 2||An||||n|| < \\Auf + \\uf, 



(u^ G{^{x-iY + {y-r^Y)u\i, v)dCd7j) 
u\[ G{^{x - if + {y- r,f)u\i, r,)didr,)dxdy 



D JD 



>b u {x,y)dxdy / u {i,r])didr] = b\\u\\ . 

JD JD 



Then from ([T^) we get 



j^Wuf + 2(a - l)||nf + 26||Mf < 0. (11) 



It is easy to see that if a > 1, i.e., /U < 0, then all solutions approach zero in L . We will 
not consider this simple dynamical case. In the rest of this paper we assume that /x > 0, 
i.e., Q < 1. 

Thus we have, for any constant e > 0, 



dt 
or 



^ "uf + 2e\\uf + 2(a - 1 - e)\\uf + 2b\\uf < 0, (12) 
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So 



By the usual Gronwall inequality (||l^) we obtain 



d ,, 1,9 ,, 1,9 (a — 1 — e)^ , 
^M^ + 2.M^< ' (14) 



\\uf < \\uofexp{-2et) + ^'^ . (15) 

When e = l — a = /i, we get the optimal or tight estimate 

\\uf < \\uofexp{-2nt) + ^. (16) 

This completes the proof of Lemma |l|. ■ 
Moreover, higher order derivatives of u are also uniformly bounded. 

Lemma 2 Suppose u is a solution of Then ||Vu|| and ||An|| are uniformly (in 

time) bounded. 

In order to prove this lemma, we recall a few useful inequalities. 
Uniform Gronwall inequality ([17|). Let g,h,y be three positive locally integrable 
functions on [to,+oo) satisfying the inequalities 

Tt-'y^^' 

with fl'^'^ gds < oi, //^"^ hds < 02 and Jl~^^ yds < 03 for t > to, where the ai(i=l,2,3) are 
positive constants. Then 

y{t + 1) < {0.2 + a3)exp{ai),foT t > to. 



Gagliardo-Nirenberg inequality (||Tl|). Let w £ L'^ CiW"^''^ (D) , where 1 < q,r < 00. For 
any integer j, < j < m, ^ < A < 1. 

\\D^w\\p < Co\\w\\l-^\\U'^w\\i: 

provided 

1 j ^ , 1 m , 1 - A 

- = - + A + , 

p n r n q 

and rn — j — J is not a nonnegative integer If m— j — ^ is a nonnegative integer, then the 
inequality (^) holds for A = ^. 

Poincare inequality (Q). For w G Hq{D), 

h\\wf < ||Vw|p, 
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where Ai is the first eigenvalue of —A on the domain D, with zero Dirichlet boundary 
condition on dD. 

Proof of Lemma |^. Due to the boundary condition @ on and the Poincare in- 
equality, we get ||Vti|p < A]^"'^||Au|p. Hence it is sufficient to prove that ||Ati|| is bounded. 
We first show that J^*^^ ||An|p(is is bounded. In fact, using 

2|(Ati,u)| < 2||An||||n|| < ^||Au|p + 2\\u\\^ , 

in (|lO|), we get 

\uf + IIAiif + 2(a - 2)||Mf + 26||uf < 0. (17) 



Since 



dr 



2b\\u\\^ + 2(a - 2)\\uf = h\\uf + 2/3(||nf + \ ^ ||n|| 

,„ „2 „2 2a -4-/3,2 (2a-4-/?)2 

= h\\uf + 2f3{\\uf + ^2 ^ 



> b\\u 



4/3 ' 
(2a -4-/3)2 



we conclude 



+ ||A„|P + < P"-"-^'' = Vl±lt±f!l^ (18) 

dt" " " " " " - 8/3 8/3 ^ ' 



Integrating (|18| ) with respect to t from t to t + 1 and noting Lemma Pj, we see that 



JI'^^ IIAnpds is bounded. 



Now, multiplying (^ by A^u and integrating over it follows that 

IIAnlP + IIA^nlP + 2 / AnA^udxd-y + all AnlP 

2dt Jd 

+ J (uj G{^{x - + {y- T]Y)u^{i,7])did'q)A^udxdy = 0. (19) 



Note that 



and 



2| / AnA\(ix(iy| < -IIA^uf + 2||Anf , (20) 
7d 2 



{u G{^{x - + {y- 7] f)u\^, r])d^drj)A\dxdy\ 
^(A«)2(^ G(^(x - + (y - ^?)u\i,r^)didrj)dxdy 
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+ ^ uAuiJ^iAGi^ix - + iy- vy)))u\t v)d(dr])dxdy 
+2 j^VuAu{jjyG{^{x - + {y- r])-^)))u\C,i])d(dr])dxdy\ 
< {a\\uf + 2A^^||VG||oohf + ^l|AG||oo||nf )||An||2 + i||AG||oo||^xf 

< (a + 2xf^Ki + ^K2)\\uf\\Auf + ^K2\\uf, 

where a,Ki,K2 are various upper bounds of G defined in (^), and R is the bound of 
the solution u as in Lemma 1. Hence by ([T9|) we get 

-^IIAnf < 2[{a + 2X~~^ Ki + ^^2)11^^11' - a + 2]\\Auf + K2\\ut ■ (21) 

(Xv Zi 

Finally, applying the uniform Gronwall inequality (pT| ) and noting Lemma |l], we conclude 
that ||Ati|p is uniformly bounded for all t > 0. This proves Lemma ||. ■ 
We now have the following global existence and uniqueness result. 

Theorem 1 Let UQ{x,y) £ L^{D) and G satisfies then the initial- boundary value 
problem (§),(!), (|) has a unique global solution u G L°^{0,oo] Hq{D)). Moreover, the 
corresponding solution semigroup S{t), defined by 

u = S{t)uo, 

has a bounded absorbing set 

Bo = {ue Hl{D) : {\\uf + \\Vuf + \\Auf)l < i?}, 

where R is a postive constant which depending on the uniform bound of \\u\\, ||Vn||, ||An||. 
Finally, the solution semigroup S{t), when restricted on Hq{D), is continuous from Hq{D) 
into Hl{D) fort>0. 

Proof. The global existence, uniqueness and absorbing property follow from standard 
arguments (e.g., |jl^) together with Lemmas ||, ^ above. The absorbing property also 
follows from these two lemmas. 

We now prove that S{t) is continuous in H"^ {D)r\HQ{D) . Suppose that uq, vq G Hq{D) 
with ||Ano||, HAvoH < 2i?i, we denote by u{t),v{t) the corresponding solutions, i.e., u{t) = 
S{t)uQ,v{t) = S{t)vQ. Let w{t) = u{t) — v{t). Then w{t) satisfies 

wt + A'^w + 2Aw + aw + w J G{^J (x - ^)2 + {y - r/)2)n^(^, r])d£,dr]+ 
V G(^(x-e)2 + (y-r?)2)(n(e, r,) + v{^, r,))w{^, v)d^dr, = 0. (22) 
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Applying the Gagliardo-Nirenberg inequality 



||^t||oo < C'o||Aii||, 

and the Poincare inequality 

„ „ 1 „ A „ 

\\w\\ < —\\Aw\\, 
we obtain (similar to the proof of Lemma ^ , 

which implies that ||Ait;(t)|p < || Awolpexp(Cit) for some positive constant Ci. This 
shows that S{t) is continuous. ■ 
This theorem implies that @-(^ defines an infinite dimensional nonlocal dynamical 
system. 

In the rest of this section, we consider the global attractor for the nonlocal dynamical 
system (Q)-(^). We will establish the following result about the global attractor. 

Theorem 2 There exists a global attractor A for the nonlocal dynamical system @, (|5|), @ . 
The global attractor is the uj— limit set of the absorbing set Bq (as in Theorem^, and it 
has the following properties: 

(i) A is compact and S{t)A = A, for t > 0; 

(ii) for every bounded set B C Hl{D), lim d{S{t)B,A) = 0; 

(iii) A is connected in Hq{D), where d{X,Y) = sup inf ||x — vWh^ii)) is the Hausdorff 
distance. 

Moreover, the global attractor A has finite Hausdorff dimension dniA) < m, where 



m~C(l + ^/i+(2a-6)^), 

where C > is a constant depending only on the domain D, and a > 0, 6 > are the 
upper, lower bounds of the kernel G, respectively. 



Proof. The existence and properties of A are quite standard now (see |17| and refer- 
ences therein). We omit this part, and only estimate the dimensions below. 



As in 1 17 1, we may use the so-called Constantin-Foias-Temam trace formula (which 
works for the semiflow S{t) here) to estimate the sum of the global Lyapunov exponents 
of A. The sum of these Lyapunov exponents can then be used to estimate the upper 
bounds of ^'s Hausdorff dimension, duiA). To this end, we linearize equation (Q) about a 
solution u{t) in the global attractor to obtain an equation for v{t) and then use the trace 
formula to estimate the sum of the global Lyapunov exponents. Doing so, we obtain 

vt + L{u{t))v = (23) 
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where 



L{u{t))v = A^v + 2Av + av + v J G{^J {x - + {y - v)dCd7] 



+2u {x - if + {y- i]f)u{i, i])v{^, r])d^dr]. 



This equation is supplemented with y, 0) =^(x,y) G Hq{D). Denote by y), . . . , ^m,(2^, y), 
m hnearly independent functions in Hq{D), and Vi{x,y,t) the solution of (|23|) satisfying 
Vi{x, y, 0) = (,i{x, y), i = 1, . . . ,m. Let Qmif) represent the orthogonal projection of Hq{D) 
onto the subspace spanned by {vi{x, y,t), . . . , Vm{x, y, t)}. 

We need to estimate the lower bound of Tr{L{u{t)Qm{t))), which gives bounds on the 
sum of global Lyapunov exponents. Note that in |jl^, the linearized equation like ( ^3|) 
is written as vt = L{u{t))v and in that case one needs to estimate the upper bound of 
Tr{L{u{t)Qm{t))). Suppose that (j)i{t), (j)mit) is an orthonormal basis (\\(pj\\ = 1) of the 
subspace Qm{t)HQ{D) for any t > 0. 

Now we estimate the lower bound of Tr{L{u{t)Qm{t))) . It is easy to see that 

Tr{Liu{t)Qm)) 



m 

= ^(A2</,,+2Ac/>, / Gi^{x-0^ + iy-vrWii,v)dCdvAj) 

3=1 ■'^ 

m 

+ Y.{2u / G{^{x-if + {y-iif)u{i,r^)<Pjdid,]A3)- 

j=l •'D 

Since (2A</)j,0j) > — (i||A0j|p + e||(/)j|p) for any constant e > 1, we get 

m -1 

Tr{L{u{t)Q^)) > Y^lil - -)\\A<t>,f + b\\4>jf\\uf + (q - e)U,f] 

3=1 ^ 

+ f22 f u(l)jdxdy{[ G{^ {x - 0^ + {y- v?Ht r])cpjdidii) 



3 

m 



-I III, 

> ^(1 - -)||A0,f + E(6||n||2 + a-e- 2a\\uf) 
j=i ^ i=i 



m 1 

: ^(1 - -)||A0,||2 + [l-ij-e+{b- 2a)\\uf]m. (24) 
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We introduce notation f{x,y) = J2 Note that m = Jj^ f{x,y)dxdy. By the 

generahzed Sobolev-Lieb-Thirring inequality ([17|, page 462), 



f{x,y)dxdy<KoY.\\AcPjf, 

where Kq > depending only on the domain D. Moreover, due to the fact that L^{D) 
L\D), 



= (/ f{x,y)dxdyY < / r{x,y)dxdy 
Jd Jd 



<KoC2j2\\A<Pjf 

m 

= Cj2\\Acl>,f 
i=i 

for some constants C2 > 0, C > depending only on the domain D. 
Thus 



1 1 1 

(l--)^||A,^,f>(l--) (25) 



Therefore, by (p3)-(pl) we have 



1 - i 

Tr{L{u{t)Q^)) > —-^m^-{fi-l + e + {2a-b)\\u\\'^)m 

> ^m^ - (/i - 1 + e + (2a - 5)^)m 

> (26) 

whenever 



m> ./[//-! + € + (2a -6)^]^. (27) 



The right hand side of (p^) has the minimal value of 



m~C(l + ^/;U + (2a-6)^) (28) 



when e = 1 + + (2a - 6)f . 

As in , we conclude that the Hausdorff dimension of A is estimated as in (|28|) . This 
proves Theorem ^. ■ 
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3 Local Swift-Hohenberg Model 



Similarly, for the two-dimensional local Swift-Hohenberg equation (^), we can obtain the 
existence of the global attractor A. We omit this part and will only estimate the dimension 
of A 

Theorem 3 There exists the global attractor A for the local dynamical system 

The Hausdorff dimension of A is finite, and dniA) < mi ~ C(l + ^/Ji), where C is 
a constant depending only on the domain D. 

Proof. As in the proof of Ttheorem we consider the linearized equation of @, 
defined by 

vt + Li{u{t))v = Q, 

where 

Li{u{t))v = A'^v + 2Av + av + Sn^v. 

Then we estimate 

Tr{Li{u{t)Qm)) 

m 

= Xl(AVj + 2A(/.j + acpj + 2,u^(l)j,(j)j) 

m 

m 1 m 

>Y^{1 - -)\\Act>,f +Y^{a - el 

j=l ^ 3=1 



m 



where we have used the fact that 3{u^ , cfij) > 0. Noting again that < C X] ll^i^jp 
and a = 1 — ^, we have 



1 - i 

Tr{Li{u{t)Qm)) > —^m^ - - 1 + e)m > (29) 



whenever 



m> ^(/i-l + e)^. (30) 

The right hand side of (^) has the minimal value of 

m ~ C7(l + ^) (31) 
when e = 1 + This completes the proof. ■ 
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4 Discussions 



In this paper, we have discussed the Hausdorff dimension estimates for the global attractors 
of the two-dimensional nonlocal and local Swift-Hohenberg model for Rayleigh-Benard 
convection. 

The Hausdorff dimension for the global attractor of the nonlocal model is estimated 

as 

m ~ C(l + ^fi + {2a-b)i^), 
while for the local model this estimate is 

m~C7(l + ^), 

where C > is an absolute constant depending only on the fluid convection domain, and 
fi > measures the difference of the Rayleigh number from its critical convection onset 
value. Note that a,b > are the upper and lower bounds, respectively, of the kernel G of 
the nonlocal nonlinearity in @. 

The two dimension estimates above differ by an absolute constant (2a — b)^, which 
depends only on the the Rayleigh number through /u, and upper and lower bounds of the 
kernel G of the nonlocal nonlinearity. Moreover, if the kernel G is a constant function 
(thus, a = b = G), then the dimension estimate for the nonlocal model becomes 

m ~ C(l + ^/2il), 

which still differs from the dimension estimate for the local model by a constant depending 
on the Rayleigh number through /x. 
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